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Abstract We report a method of solving for canonical scalar field exact solution in a non-flat FLRW 
^ ' universe with barotropic fluid using non-linear Schrodinger (NLS)-type formulation in comparison to 

O^' the method in the standard Friedmann framework. We consider phantom and non-phantom scalar 

j ' field cases with exponential and power-law accelerating expansion. Analysis on effective equation of 

bi)! state to both cases of expansion is also performed. We speculate and comment on some advantage and 

disadvantage of using the NLS formulation in solving for the exact solution. 
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1 Introduction 

I In the past decade, it has been observed that universe is now in acceleration phase [Tl[51[31|3] while 

inflationary scenario of the early universe [5] is strongly confirmed by cosmic microwave background 
data [HITIIHIIS]. In both circumstances, the universe experiences accelerating expansion which can be 
, attained by exploiting some dynamical scalar field with time-dependent equation of state coefficient 

^ • w^{t) < —1/3, or a cosmological constant with wa — —'^ [ID]- It was further suggested that the 

I scalar field could have phantom behavior with < — 1 as one considers kinetic energy term in its 

Lagrangian density to be negative [TT|fT2] . Strong supports to the phantom idea are from observations 
previously, e.g. combined cosmic microwave background, large scale structure survey and supernovae 
type la without assuming flat universe yields w^^ = — l.OB^to og 15] while using supernovae data alone 
assuming flat universe, w^^o = — 1.07±0.09 [13j. The subscript denotes the value at present. Moreover, 
most recent WMAP five-year result [7l[8] combined with Baryon Acoustic Oscillation (BAO) of large 
scale structure survey from SDSS and 2dFGRS [14] and type la supernovae data from HST [3\ , SNLS 
[3] and ESSENCE 13. assuming dynamical w with flat universe yields —1.38 < w^^ < —0.86 at 95% 
confident level and w^^ = —1.12 ±0.13 at 68% confident level 0. With additional BBN constraint of 
limit of expansion rate [T5lll6j. —1.32 < w^^o < —0.86 at 95% confident level and w^.o = —1.09 ± 0.12 

This work is supported by a TRF-CHE Research Career Development Grant of the Thailand Research Fund 
and the Naresuan Faculty of Science Research Scheme. 

Burin Gumjudpai 

Fundamental Physics & Cosmology Research Unit, The Tah Foe Academia Institute (TPTP) 

Department of Physics, Naresuan University, Phitsanulok 65000, Thailand 

and 

School of Physics, Institute of Science, Suranaree University of Technology 
111 University Avenue, Nakhon Ratchasima 30000, Thailand 
E-mail: buring@nu.ac.th 



2 



at 68% confident level. This suggests that phantom field has firmed status in cosmology. However, 
phantom field does result in unwanted Big Rip singularity in a flat FLRW universe [17] . This raises up 
many attempts to avoid the singularity based on both phenomenological and fundamental inspirations 

Recently, there are a few proposals for mathematical alternatives to the conventional Friedmann 
formulation of canonical scalar field cosmology, such as non-linear Ermakov-Pinney equation [l9l[20] . 
Moreover, a non-Ermakov-Pinney equation for the same system was also proposed in form of a non- 
linear Schrodinger-type equation (hereafter-NLS) and it was found that solutions of the NLS-type 
equation are correspondent to solutions of the generalized Ermakov-Pinney equation [20..21p ]. Conclu- 
sion of how to relate NLS quantities to quantities of standard Friedmann formulation is shown in [23] 
which gives extension to phantom field case. It also shows that the NLS wave function is in general non- 
normalizable. Expressing cosmological quantities in form of NLS quantities may suggest an alternative 
way of solving problems in scalar field cosmology. In such method, presumed knowledge of scale factor 
function with time a(t) must be given first and later one can evaluate NLS potential based on a(t) 
assumed. Here we attempt to solve for scalar field exact solution within the NLS framework in various 
cases with flat and non-flat spatial geometries. We compare the results to the solutions obtained in 
standard Friedmann formulation. Exponential expansion a ^ exp(t/r) and power-law expansion a 
are assumed where r are finite characteristic time and g is a positive constant. 

This article is organized as follow. We express our cosmological system in Sec. [5] before introducing 
the NLS formulation in Sec.[3l Afterward, we consider each model of expansion separately. Exponential 
expansion is presented in Sec. [4] where we obtain exact solution from effective equation of state and later 
we solve for exact solution from Friedmann system. Afterward, in Sec. [5l we consider NLS formulation 
for the exponential expansion and solve for exact solutions in NLS framework. We analyze effective 
equation of state for the exponential expansion case in Sec.|6l When considering power-law expansion, 
we work and organize the contents in the same spirit and order as in the previous sections. Beginning 
from Sec. [T] [8] and m Finally, conclusion is made in Sec. [TOl 



2 Cosmological System 

To be realistic, two perfect fiuids are considered in our system: barotropic fiuid and scalar field fiuid. 
The perfect barotropic fiuid pressure p~f and density obey an equation of state, — {'^—l)p-y — w^pj 
while for scalar field, = w^p^. Total density and total pressure are ptot = P-( + P4> and ptot = P-y +P<t>- 
The effective equation of state is weighed- value of these two components, 

p^w^ + p^w^ 

WcS = . [1) 

Ptot 

For the barotropic fiuid, its equation of state coefficient w-y is written in term of n. We set w-y = (n— 3)/3 
so that n = 3(1 + w-y) — 37, hence Wy — —1 corresponds to n = 0, Wy — —1/3 to n = 2, w-y = to 
n = 3, Wy ~ 1/3 to n — 4, and Wy — I to n — 6. The conservation equation is therefore 

p~f = —nHpy (2) 

with solution obtained directly, 

P7 = ^ . (3) 

therefore p^ = [{n — 3)/3](£'/a") , where a is scale factor, the dot denotes time derivative, H = d/a is 
Hubble parameter and _D > is a proportional constant. The scalar field considered here is minimally 
coupling to gravity with Lagrangian density, C — {l/2)e(jp — V{(j)) , where e = 1 for non-phantom case 
and —1 for phantom case. Density and pressure of the field are given as 

P0 = + V{(^) , P4, = ie02 _ y(^) ^ (4) 



^ Considering Bianchi I scalar field cosmolog y, on e can also construct a corresponding linear Schrodinger-type 
equation by redefining cosmological quantities |22) 
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therefore 



P<f> 
= — 



e(/)2 - 2F((/>) 



P4> e</)2 + 2F(0) 



The field obeys conservation equation 



dV 



= 0. 



(5) 



(6) 



Considering Friedmann-Lemaitre-Robertson- Walker (FLRW) universe, the Friedmann equation and 
acceleration equation are 



H — -z-ptot K 



a 
a 
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■Ptot(l + iWcs) 



(7) 

(8) 



where = 8ttG = 1/Mp, G is Newton's gravitational constant, Afp is reduced Planck mass, k is 
spatial curvature. Using Eqs. Q, ([U, ([6]) and (O, it is straightforward to show that 



nD 

2k 
3^ 



n-Q\ D 
6 / a" 



(9) 

(10) 



3 Non-linear Schrodinger-type formulation 

Correspondence between non-linear Schrodinger formulation for canonical scalar field cosmology with 
barotropic fluid was shown in [21j and was concluded recently in |,23j . In the Schrodinger formulation, 
wave function u{x) is related to scale factor in cosmology as 

EE a(t)-"/2 , (11) 

and Schrodinger total energy E and Schrodinger potential P{x) are linked to cosmology as 

^---TT^' (12) 

P{x)^'-^a{tT4{tf. (13) 
These quantities satisfy non-linear Schrodinger-type equation: 

^n{x) + [E- P{x)\ u{x) = -!^u(x)(4-»)/" . (14) 

The mapping from i to a; is via x — u{t), such that 

x{t)^u{x), (15) 
0(i) = ^{x) . (16) 

We comment that the relation t/;'{x)'^ ~ {^/nK^)P{x) in Ref [21] does not include phantom field case. 
Modification is made in recent work [23 so that the solution includes the phantom field case, therefore 



^(^)^± 2 /• Wdx. (17) 



If P{x) and n ^ 0. There exists an inverse function of ip{x) as ip ^{x). Therefore x{t) = ip ^ o (j)(t) 
and the scalar field potential, V o a~^{x) can be expressed as function of time, 

, 12 fduV 2u2 12^2 3fcu4/n 
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4 Exponential expansion 

4.1 Solution solved from effective equation of state for A: = case 
Exponential expansion reads 

a{t) = exp {t/r) , (19) 

where r is a positive constant. Flat universe undergoes exponential expansion only when Wcs = — 1. 
The effective equation of state, (Eq. ([T])) with Eqs. (jl]) and ^ can therefore be written afl 

e0' = -^P7: (20) 
which can be integrated directly, using Eq. to 

(/.(t) = ±2rA/^e-"*/2- + 0o. (21) 
V on 

The solution above is obtained when assuming phantom scalar field, i.e. e = — 1. If the scalar field is 
not phantom, the solution is imaginary. 



4.2 Solution solved from Friedmann formulation 



Another way to find the exact solution is to use Eq. (fT51) . in Eq. Therefore 

2k nD 



6,^(^)^ = 9e-2*/--^e""*/^ (22) 



which gives an integration: 




0(<)=± / jif^e-2*/--^e-W^)dt. (23) 



4-2.1 Simplest case 



In the case of fc = and D — 0, the integration yields a constant (f)Q. Eq. ^ becomes = —1. This 
is a cosmological constant as seen in simplest model of exponential expansion. When assuming only 
k = and e = — 1 but with D ^ 0, the solution of Eq. is the same as the Eq. (pi]) previously. For 
a scalar field domination in a non-flat universe (_D = 0, fc 7^ 0), the solution is 

=±I./2^e-t/- + 0„. (24) 
K V e 

where k and e must have the same sign, otherwise the solution is imaginary. 



4-2.2 The case of non-zero k and non-zero D 



When k and D are both not negligible. Performing integration to the Eq. ([23]) is more complicated 
and could be impossible unless assumption of barotropic fluid type. When assuming a particular type 
of fluid in the integration, i.e. n = 0,2,3,4 and 6, analytical solution can be found for all n vales in 
complicated forms. For example, the simplest among these is dust case (n = 3) which has solution: 

with additional rule that fc > and e = 1 otherwise it is imaginary. In the next section, we will show 
how to obtain solution in NLS formulation for dust and radiation cases. 

^ We are not considering a cosmological constant but a dynamical scalar field and a barotropic fluid which 
together yield WeS = — 1- 
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5 Exponential expansion: Solutions solved with NLS formulation 

For exponential expansion, following Eqs. and (|15p . we get 

u{x) = x{t) = exp (— nt/2r) . 
Integrating the above equation, hence parameters x and t scale as 

2r 



where is an integration constant. The reverse is 

2t 



t{x) 



In [(-n/2T) (a; - xo)] , 



where the condition x < xq must be imposed. Now we can write wave function as 



u{x) 



2t 



{x - Xo) 



which is a linear function. Using Eq. ([2^ . hence the Eq. (fT5|) reads 



Here the Schrodinger kinetic energy term is 

m 



kn 



,(n-2)t/T 



Expressing in Schrodinger formulation, these functions are written in term of x, 

-2(«-2)/n n'^n'^D 



12 



Tlx) 



kn 

T 



2t 



(a; - xq) 



-2(n-2)/n 



(26) 
(27) 

(28) 
(29) 

(30) 

(31) 

(32) 
(33) 



In order to obtain the scalar field potential V{t), we use Eqs. (|lip . (|12p . (|13p in Eq. (|18p . we finally 
obtain 



Vit) 



2k 



-2t/T 



n — 6 
6 



which is checked by using Eq. (fT9|) in standard formula pO]) . We use Eq. ([32]) in Eq. (|17p . then 

±2 



V'(a;) 



27 



2t 



(x - xo) 



-2(n-2)/n 



n'^n'^D 
12e 



dx . 



(34) 



(35) 



We will integrate this equation in cases oi k — and k ^ 0. 



5.1 The case k^O 

When fc = and D integrating Eq. ([35|) and transforming x to t with Eq. ([27|) yields same result as 
Eq. ()2ip obtained by solving effective equation of state equation or by integrating from the Friedmann 
formulation. Real solution exists only when the scalar field is phantom. With the solution pT|) . the 
scalar field potential in term of (j), reads 



n — 6 \ 3n 



6 



4t2 



(36) 
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5.2 The case k 

When k ^ and D ^ 0, the integral ([35|) can be integrated yielding complicated hypergeometric 
function even when n is not specified. The case n = is excluded from our consideration by the reason 
mentioned in Sec[3l For naturalness, we consider radiation [n = 4) and dust {n = 3). 

5.2.1 Radiation case 

Radiation fluid corresponds to n = 4, the integral ([55]) becomes 

1 




kr 



1 



e {x- xq) 



dx . 



(37) 



Here x could be negative, e can possibly be either ±1. The solution in radiation case is 



■0(2;) 



4 kr 

--D{x ~ xqY -[x - xo) 

3 



, fcr / 3 

±-— rW -— arctan ■ 
4k2 V 



^K^D{x~XQ)l-ie\ +kT/e 



[AK^/D/{e^/i)]^~[AK^D{x - a;o)V3] - fcr(a; - xq) 



- V'o 



(38) 

allowing only e = 1 case for the solution to be real. Transforming x scale to the t scale using Eq. (P7|) . 
the solution therefore reads 



p — it/T I p — 2t/7 

3 ' 2k2 




, fcr / 3 
± — -— X arctan ■ 
4k2 V De 



iAK^DT/{3e)]e-^'/^ + kT/e 



(39) 



The solution above, when assuming fc = 0, reduces to the solution (j2ip when n = 4, confirming the 
correctness of the result obtained. 



5.2.2 Dust case 

The integral IpS]) in the dust case n — Z reads 

2 



with solutions 



ip{x) 



i^{x) = ±1 




1/3 



1 



3k^D 



3t2 
2 



1/3 



4fc 



e {x — xo)^/'^ 

(.T-Xo)^/^ 



4 e 

3/2 



dx , 



With similar procedure to the radiation case, using ((27|) . the solution is therefore, 



2t 



2fc 



3/2 



(40) 



(41) 



which is the same as Eq. ()25p derived from Friedmann formulation. This solution when assuming fc = 
is exactly the same as the solution (|2ip when n = 3 (dust fluid). This also confirms that our results 
from NLS formulation are correct. The NLS solution can solve the case when fc and D are non-zero 
together without knowing n value while the standard procedures in Sec. 14.21 can not unless assuming a 
particular value n = 0, 2, 3, 4, 6. However, it must be noticed that one can not reduce the NLS solutions 
([39|l and ([42|) to the D — case directly since there are mixed multiplication term of n and fc in the 
solution and also the value of n has already been put in. Hence setting 13 = in ((39|) and (|42l) can not 
be considered as a pure scalar field dominant case. 
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6 Exponential expansion: Analysis on effective equation of state coefficient 

The exponential expansion in our scenario is caused from mixed effect of fluids and spatial curvature. 
We discuss mixed effect on equation of state here. Definition of effective equation of state coefficient, 
WcS = {ptj^Wcp + P-yW-y)/ptat together with Eqs. ([4]), (|22p and (|34p in context of exponential expansion 
becomes 

-1 - (fcrV3)e-^*/- 
"'off = 1 + kr^,-2t/r ' (43) 



which is infinite when 



t^^H-kr'). (44) 



Infinity can possibly happen only when k = —1 because logarithm function forbids negative domain. 
In order to acquire exponential expansion in flat universe, one needs to have WcS = — 1, but this is not 
true when k term is non-trivial. Therefore we can only express Wrf, in term of Wcs as 

^ [(3fc^^)e-^*/" + 3/(A.V)]u;eff ~ [(n - 3)/3] j^e'"^/" 
'^'^ (3/fc/K2)e-2t/r + 3/(K2r2) - De""*/^ ' ^ ' 

for exponential expansion. The Eq. does not depend on properties (n) or amount {D) of the 
barotropic fluid. It reduces to Wcs — —I when fc = as expected. Considering Eq. ([^5]) . if D = and 
fc = 0, it yields = Wcs while setting D = alone also gives the same result. 



7 Power-law expansion 

7.1 Bound value of (t){t) from effective equation of state for fc = case 
In power-law expanding universe, scale factor evolves with time as 

ait) = . (46) 

where q > is a constant. In flat (fc — 0) universe, it is known that the power-law expansion, is 
attained when —1 < WcS < —1/3 where q — 2/[3(l -I- Wcs)]- The effective equation of state, (Eq. ([T|)) 
hence is a condition 

i.e. < eij)'^ + < (2/3)ptot- Both values of e can be assigned and the power-law expansion is 

sustained as long as the condition is satisfied. 

7.2 Solution solved from Friedmann formulation 

If we directly consider Eq. ([9|) , the solution for power-law expansion is an integration: 



^(,)=± / Jif4^ + 4^_^)dt. (48) 

' ^' e I k2<2 K^t^q 3^9" ' ^ ^ 
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1.2.1 Simplest case 

Simplest integration case is when fc = and D = 0. The solution of Eq. (f48|) is well known [24] 

= ±\/^lni + 0o, 
V eK^ 



(49) 



provided that q and e have the same sign. Considering power-law inflation, the WMAP five-year 
combined analysis based on flat and scalar field domination assumption yields g > 60 at more than 99 
% of confident level otherwise excluded while q ^ 120 is at boundary of 68% confident level [9]. These 
results base on single field model which we can applied the above solution to. When assuming only 
fc = with -D 7^ 0, the solution of Eq. (gHl) is 




{l + y/l - {nDn'^lQq) t-9"+2^ ' 




l-qn+2 _^ 



uDhl 



where, when q = 2/n, the field has infinite value. The last logarithmic term in the bracket is an 
integrating constant which is valid only when q < 0. To attain power-law expansion, q must be 
positive. Hence, this term is not defined for power-law expansion. We will see later that the NLS result 
does not have this problem. For the reverse case, £> = 0, fc 7^ 0, the solution is 






t-2q+2 + 1 I + ,^0 , (51) 



which becomes infinite when q = 1. The values of q, k and e must have the same sign in all terms of 
the solution otherwise becoming imaginary. Hence, for q > 0, the condition for the solution to be valid 
is fc 1 and e — 1. 

7.2.2 The case of non-zero k and non-zero D 



(50) 



When considering non negligible value of both fc and D, the Eq. can not be integrated analytically 
except when setting n = 2 (w^ — ^1/3) which is not natural fiuid. Hence it is not considered. 



8 Pov^rer-law expansion: Solutions solved with NLS formulation 



Power-law expansion cosmology in NLS-type formulation is presented and concluded in |23j . Important 
functions needed for evaluating the field exact solutions are 



X = a(t) = — +Xq, 



t{x) 

P{x) 
Vit) 



1 



[-P{x~xo)f^ ' 
2q 2k uD 



kH"^ 3t9" 

2qn 1 
(qn — 2)2 (x — xq)^ 
q{3q~l) 2k 



2^2 



^2f-2q 



We use Eq. ^ in Eq. ([T7D, then 
±2 



kn 
~2 



-2 



{qn -2){x- xq) 
n-6\ D 
6 / ■ 



2q{n-2) j (qn-2) 



12 



■0(a;) 



2qn 



e{qn — 2)^ (x — Xq) 



kn 
2+27 



-|29(n-2)/(g«-2) 



(qn -2) {x- xo) 



I2e 



(52) 
(53) 
(54) 
(55) 
(56) 

Ax .(57) 
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We consider the solution in cases of fc = and k ^ 0. Recall that setting D — can not be considered 
as an absence of barotropic fluid due to existence of n in the other terms. 



8.1 The case k = Q 

Solution to the integral ([57|) for fc = case is 



■)p{x) = ±. 



eK'^{qn — 2)^ 



K?Dn{qn — 2)^ 
24^ 



{x - 



In 



1 + ^1 - [ti'^Dniqn - 2)2/24(7] {x - Xq)'^ 4gn 



{x ~ xo) 



e{gn - 2f 



Transforming to the t variable using Eq. (j52p . we obtain, 



m = ± 



1 / 2q 



In 



[l + v^l - (nL>K2/6q)i-9»+2 



2^|1- I — ^ I <-9"+2 + In ' ^ 



V 6g 



2qn 



This solution differs from the solution ([50]) only the last logarithmic term in the bracket which is only 
an integrating constant. When q — 2/n or n = 0, the field has infinite value. The last logarithmic term 
does not restrict the sign of q. Only q and e must have the same sign for the solution to be real. 



8.2 The case k^O 

In case of non-zero k and non-zero D, the integral ([57| can not be integrated analytically even when 
assuming each n value except when n = 2 which is not natural fluid. 



9 Power-law expansion: Analysis on effective equation of state coefficient 

Similar to the analysis in Sec. [SI mixed effect of the two fluids and spatial curvature results in power-law 
expansion. The coefficient WeS, with Eqs. ([3]), ([5^ and ([55]) . reads 

(-3g2 + 2g)t2^-2^fc 

^'eff = 3g2,2,-.+3fc ' (60) 



1/(29-2) 



which becomes infinity if 

We can also express in term of Wos as 

^ [(3gV«2)^-2 ^ (^sk/K^)t-^i]w,« - [jn - 3)/3] Dt-'J'' 

for power-law expansion. The Eq. (|62p . when _D — and fc = 0, yields — WcR as expected. Similar 
to the case of exponential expansion, setting D — Q alone also yields — Wcs- In flat universe, 
power-law expansion happens when Wcs lies in an interval (—1, —1/3). But in fc 7^ universe, it is no 
longer true. Considering flat universe, setting fc = in Eq. ([60| yields q — 2/[3(l + ti;off)]. The condition 
— 1 < Wcff < — 1/3 therefore corresponds to g > as known. The condition also yields 

+ <w^ < -r-(\ + w^\P^. (63) 

P<h 3 V3 Pa 
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If there is more non- negligible radiation fluid (with w-^ — 1/3), it is noticed that the interval shifts to 
the more left. For example, setting = O.lp^, the interval shifts to about —1.133 < < —0.4. If 
we assume more realistic situation when dust (dark matter and other matter elements) is presented. 
The dust density and dark energy is about 28% and 72% of total density, therefore p-y ~ (28/72)p^ ~ 
O.389p0, the interval is —1.389 < < —0.463 which covers valid range of recent observational data, 
assuming dynamical w with flat universe, —1.38 < w^^o < —0.86 at 95% confident level [9]. 



10 Conclusions 

This letter reports and demonstrates a method of solving for canonical scalar field exact solution in 
a non-flat FLRW universe with barotropic fluid using non-linear Schrodinger (NLS)-type formulation 
in comparison with the method in the standard Friedmann framework. We consider phantom and 
non-phantom scalar field cases with exponential and power-law accelerating expansion. We evaluate 
all NLS quantities needed to find the solution, e.g. non-normalizable wave function and Schrodinger 
potential. Our process is reverse to a problem solving in quantum mechanics that the wave function is 
expressed first by the expansion function, a{t) before evaluating the Schrodinger potential based on a 
known expansion function. In NLS formulation the total energy E is negative. We do an analysis on 
effective equation of state to both cases of expansion. We expresses Wcs in term of q and k. In a flat 
universe, in order to have power-law expansion, the interval (—1, —1/3) of the Wfj,, is shifted leftward 
to more negative if more barotropic fluid density is presented. 

Within framework of the standard Friedmann formulation, we obtained exact solution in various 
cases. Later we solved the problem using NLS formulation, in which the wave function is equivalent 
to the scalar field exact solution. NLS method is restricted by the fact that its scalar field solution is 
valid only when the barotropic fluid density is presented. Setting D = does not imply the absence 
of barotropic fluid because the barotropic fluid parameter n still appears in the other terms of the 
Schrodinger potential. Therefore NLS formulation can not be applied to situation when the scalar field 
is dominant and D ^ 0. Hence it is more suitable for a system of dark energy and dust dark matter 
fluid. This is a disadvantage point of the NLS formulation. Transforming from standard Friedmann 
formulation to NLS formulation makes n appear in all terms of the integrand and also changes fluid 
density term D from time-dependent term to a constant E. Hence the number of a;(or equivalently 
t)-dependent terms is reduced by one. This is a good aspect of the NLS. In both Friedmann formulation 
and NLS formulation, the solutions when k ^ and D ^ are difficult or might be impossible to solve 
unless assuming values of q and n. Hence reduction number of x-dependent term helps simplifying the 
integration. There are also other good aspects of NLS formulation. Firstly, in the case of exponential 
expansion with NLS formulation, the solution when fc 7^ and D ^ Q can be obtained without assuming 
n value while n — 0,2,3,4,6 must be given if working within Friedmann formulation. Secondly, for 
power-law expansion with fc = 0, the result ()59|1 obtained from NLS formulation has integrating 
constant that does not restrict q value while (|50p obtained from Friedmann formulation needs q < Q 
which violates power-law expansion condition (g > 0). For power-law expansion, the most difficult case 
is when k ^ Q with Z? 7^ 0. In both formulations, the integral can not be integrated unless assuming 
n ~ 2 (equivalent to = —1/3) which is not a physical fluid. We introduce here alternative method 
to obtain scalar field exact solution with advantage over and disadvantage to standard Friedmann 
formulation. The NLS formulation could render more interesting techniques for scalar field cosmology. 
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